IFIC/09-09 



O 
O 

(N 



LHC phenomenology of the /iz^SSM 
A. BartQ 

Institut fiir Theoretische Physik, Universitat Wien, 
A~1090 Vienna, Austria 

M. HirsctQ and A. Vicent^S 

AHEP Group, Institut de Fisica Corpuscular - C.S.I.C. & Universttat de Valencia 
Edificto Institutos de Paterna, Apt 22085, E~46071 Valencia, Spain 

S. LiebleJl and W. Porotfl 

Institut fiir Theoretische Physik und Astronomic, Universitat Wiirzburg 
Am Hubland, D-97074 Wiirzburg, Germany 



The /ii/SSM has been proposed to solve simultaneously the /i-problem of the MSSM and explain 
current neutrino data. The model breaks lepton number as well as R-parity. In this paper we study 
the phenomenology of this proposal concentrating on neutrino masses and the decay of the lightest 
supersymmetric particle (LSP). At first we investigate in detail the /ii/SSM with one generation of 
singlets, which can explain all neutrino data, once 1-loop corrections are taken into account. Then 
we study variations of the model with more singlets, which can generate all neutrino masses and 
, mixings at tree-level. We calculate the decay properties of the lightest supersymmetric particle, 

^ l' assumed to be the lightest neutralino, taking into account all possible final states. The parameter 

I ' regions where the LSP decays within the LHC detectors but with a length large enough to show 

a displaced vertex are identified. Decay branching ratios of certain final states show characteristic 
correlations with the measured neutrino angles, allowing to test the model at the LHC. Finally we 
briefly discuss possible signatures, which allow to distinguish between different R-parity breaking 
models. 

m 
> 

V/-^ . I. INTRODUCTION 

cn 

ff^ ■ The Minimal Supersymmetric extension of the Standard Model (MSSM) ilj] assumes that R-parity is 

I conserved. R-parity (Rp) 0, defined as Rp — (_i)3B-i-l-i-2S^ originally introduced to guarantee the 

On . stability of the proton in supersymmetric models 0, 0]- It has two immediate consequences: First, the 

' lightest supersymmetric particle (LSP) is stable. For cosmological reasons a stable LSP has to be electrically 

neutral, thus leading to the "standard" missing momentum signature of SUSY. Second, the MSSM with Rp, 
1 ^ ■ for the same reasons as the SM, predicts zero neutrino masses. 

r> I Neutrino oscillation experiments have demonstrated that at least two neutrinos have non-zero mass 

0, 0]. Especially remarkable is that data from both atmospheric neutrino Q and from reactor neutrino 
measurements [9| now show the characteristic L/E dependence expected from oscillations, ruling out or 
seriously disfavouring other explanations of the observed neutrino deficits. It is fair to say that with the 
most recent data by the KamLAND 9], Super-K [l^ and MINOS collaborations [llj neutrino physics has 
finally entered the precision era. (For the latest evaluation of allowed neutrino parameter regions, see for 
example the updated fits in [l2[.) 
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Non-zero neutrino masses can be easily included into the standard model by simply adding right-handed 
neutrinos, postulating the existence of a {AL — 2) dimension-5 operator Il3l of unspecified origin or by 
introducing the seesaw mechanism with either fermionic singlets [3, [iH Il6| . a scalar triplet [13, [3] or 
fermionic triplets [l^. Neutrino masses could be induced also at 1-loop-level [1^ or even at 2-loop order 

While all of the neutrino mass models mentioned above can be easily supersymmetrized, there is also an 
entirely supersymmetric possibility to generate Majorana neutrino masses: R-parity violation [23l [2^. [25|. 
Different models of (lepton number violating) R-parity violation have been discussed in the literature. Within 
the MSSM particle content R-parity can be broken explicitly either by bilinear or by trilinear terms [l^l ■ The 
huge number of free parameters in the trilinear model, however, makes such a general ansatz rather arbitrary. 
Attempts to reduce the number of free parameters based on discrete symmetries have been discussed in the 
literature [H, [l^, [H, . One could also postulate that lepton number is conserved at the superpotential 
level, broken only by the vacuum expectation value (vev) of some singlet field [s^l • This is called spontaneous 
R-parity violation (s-^p). ^ Bilinear ^p(b-^p) can be understood as the low-energy limit of some s-.^pmodel, 
where the new singlet fields are all decoupled. Such a bilinear model has only six new parameters and is 
thus more predictive than the general case with all possible bilinear and trilinear couplings. 



The phenomenology of SUSY has been studied extensively in the past, for reviews see 33, 34 [. Neutrino 
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masses have been calculated with trilinear couplings 2J| and for pure bilinear models [3£ 
angles are not predicted in either schemes, but can be easily fitted to experimental data. In bilinear schemes 
the requirement to correctly explain neutrino data fixes all ^p couplings in sufficiently small intervals such 
that in some specific final states of the decays of the LSP correlations with neutrino angles appear. This has 
been shown for a (bino-dominated) neutralino LSP in [s^ l , for charged scalar LSPs in (39j and for sneutrino, 
chargino, gluino and squark LSPs in [ioj . Such a tight connection between neutrino physics and LSP decays 
is lost, however, in the general trilinear- plu s-bilinear case. (For some recent work on collider phenomenology 
in trihnear I^p , see for example [ill. l42l. lisl. [43| and references in [ssj.) 

The superpotential of the MSSM contains a mass term for the Higgs superfields, jiHdHu- For phenomeno- 
logical reasons this parameter must be of the order of the electro-weak scale. However, if there is a larger 
scale in the theory, like the grand unification scale, the natural value of ji lies at this large scale. This is, 
in short, the u-problem of the MSSM [H]. The Next-to-Minimal SSM (NMSSM) provides a solution to this 
problem [l^H^], at the cost of introducing a new singlet field. The vev of the singlet produces the /i term, 
once electro-weak symmetry is broken. (For some recent papers on the phenomenology of the NMSSM, see 
for example [is . l49l |50| and references therein.) 

The /^I'SSM 5l| proposes to use the same singlet superfield(s) which generate the /i term to also generate 
Dirac mass terms for the observed left-handed neutrinos. Lepton number in this approach is broken explicitly 
by cubic terms coupling only singlets. Rp is broken also and Majorana neutrino masses are generated once 
electro-weak symmetry is broken. Two recent papers have studied the /ij/SSM in more detail. In [s^ the 
authors analyze the parameter space of the /ij/SSM, putting special emphasis on constraints arising from 
correct electro- weak symmetry breaking, avoiding tachyonic states and Landau poles in the parameters. The 
phenomenology of the /xi/SSM has been studied also in In this paper formulas for tree-level neutrino 
masses are given and decays of a neutralino LSP to two-body (VF-lepton) final states have been calculated 

Hi. 

We note that similar proposals have been discussed in the literature. [54] studied a model in which the 
NMSSM singlet is coupled to (right-handed) singlet neutrino superfields. Effectively this leads to a model 
which is very similar to the NMSSM with explicit bilinear terms, as studied for example also in [55|. In [5^ 
the authors propose a model similar to the /i^^SSM, but with only one singlet. 

In the present paper, we study the phenomenology of the /ij/SSM, extending previous work [sil. [s^. [53| . We 
consider two different variations of the model. In its simplest form the /xj^SSM contains only one new singlet. 



^ The first model to propose s-Ijip [23l| used the left-sneutrinos to break Rp. This leads to a doublet Majoron, now ruled out 
by LEP data [H. 

^ In [s^ l it has been proposed that the trilinear parameters follow the hierarchies of the standard model Yukawa couplings. 
This is very similar to the pure bilinear model, which in the mass eigenstate basis has effective trilinear parameters given by 
products of bilinear parameters and down quark/charged lepton Yukawa couplings. 
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This version produces one neutrino mass at tree-level, while the remaining two neutrinos receive mass at 
the loop-level. This feature is very similar to bilinear R-parity breaking, although as discussed below, the 
relative importance of the various loops is different for the explicit bilinear model and the /zi^SSM. As in the 
explicit bilinear model neutrino angles restrict the allowed range of I^p parameters and correlations between 
certain ratios of decay branching ratios of the LSP and neutrino angles appear. In the second version we 
allow for n singlets. Neutrino masses can then be fitted with tree-level physics only. However, many of the 
features of the one generation model remain at least qualitatively also in the n singlet variants. LSP decays 
(for a bino or a singlino LSP) can be correlated with either the solar or atmospheric angle, thus allowing to 
construct explicit tests of the model for the LHC. In contrast to [5^ we consider all kinematically allowed 
final states. This does not only cover scenarios where two-body decays are important, but also those where 
three-body decays are dominant. In addition we show that even in the scenarios where two-body decay 
modes in singlet Higgs bosons dominante, the lifetime can be such that the LSP decays outside the detector. 

This paper is organized as follows. In the next section we outline the model, give the soft breaking terms, 
discuss the mass matrices and calculate approximate formulas for neutrino masses. We will not use the 
approximate formulas in our numerical analysis, but give them explicitly because they allow to understand 
in an easy way our numerical results qualitatively. In Section III we discuss existing constraints on the model 
space, apart from neutrino physics, and outline the properties of the "standard" points, which we will use in 
our numerical analysis. We then turn to the collider phenomenology of the model. In Section IV we study 
the one generation variant of the /ii^SSM. Decays of scalars are briefly discussed, before calculating decay 
properties of the neutralino LSP. Section V gives a discussion of the LSP phenomenology for the n generation 
variant, although we will mainly focus on two generations. Similarities and differences to the one generation 
model are discussed. In Section VI we then give a short, mostly qualitative discussion of possible signals 
which might give some hints which R-parity breaking model is indeed realized in nature, before closing with 
a short summary. Mass matrices and couplings are given in various appendices. 

II. MODEL BASICS 

In this section we introduce the model, work out its most important properties related to phenomenology 
and neutrino masses and mixings. As explained in the introduction, we will consider the n generations case 
in this section. Approximate formulas are then given for scalar masses for the one (1) ^''^-model and for 
neutrino masses for the 1 and 2 i>'^-model. 



A. Superpotential 

The model contains n generations of right-handed neutrino singlets. The superpotential can be written as 

W - h!iQ^U,Hu + h%Q,D,Hd + h%L,E,Hd 

+ hl'L,V^Hu-Xs^sHdHu + ^K,tn^^'',V^K ■ (1) 

The last three terms include the right-handed neutrino superfields, which additionally play the role of the 
$ superfield in the NMSSM |4f||, a gauge singlet with respect to the SM gauge group. The model does not 
contain any terms with dimensions of mass, providing a natural solution to the /i-problem of the MSSM. 
Please note, that as the number of right-handed neutrino superfields can be different from 3 we use the 
letters s, t and u as generation indices for the superfields and reserve the letter i, j and k as generation 
indices for the usual MSSM matter fields. 

The last two terms in ([T]) explicitly break lepton number and thus R-parity giving rise to neutrino masses. 
Note that Kstu is completely symmetric in all its indices. In contrast to other models with R-parity violation, 
this model does not need the presence of unnaturally small parameters with dimensions of mass, like in 
bilinear R-parity brea king models [s^l, and there is no Goldstone boson associated with the breaking of 
lepton number [2^ [H^, ISSj . since breaking of Rp is done explicitly. 
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For practical purposes, it is useful to write the superpotential in the basis where the right-handed neutrinos 
have a diagonal mass matrix. Since their masses are induced by the k, term in ([T]) , this is equivalent to writing 
this term including only diagonal couplings: 

n 

KstuK^tK =^ y^l^sil^lf (2) 



B. Soft terms 

The soft SUSY breaking terms of the model are 

Vsoft = Vsoft + ^soft ■ K'^) 

^soft^^ contains all the usual soft terms of the MSSM but the B^-term 

+ mjj^HrHS + m]j^HZ*HZ ~ \ [A'hB^B" + AhW^W^ + Ahg'^g'' + h.c] (4) 



+ ea6 [Ti]Q1U*Hl + riQ\b*H^^ + T'fL\E*H^^ + h.c] 
and V^^^f'^^^ includes the new terms with singlets: 

y::f' = m^^^t^r + [niL'^^i^fHt - Tiv';mHi + h.c] + [^t^^-dididi + h.c] (s) 

In these expressions the notation for the soft trilinear couplings introduced in [H, H^l is used. Note that the 
rotation made in the superpotential does not necessarily diagonalize the soft trilinear terms T^*" implying 
in general additional mixing between the right-handed sneutrinos. 

C. Scalar potential and its minimization 

Summing up the different contributions, the scalar potential considering only neutral fields reads 

V^Vd + Vf + Vsoft (6) 

with 

Vd = \{9'+9''){\K?-m'-Y^\h?Y (7) 

i—1 

3 n 

Vf = \h^i>.i>^ - XsD^^HX + \\si^sH'J' + E K'^sH^j' + K'^^f^u - A.i/„"i7° + -Ksi^ri' , (8) 

where summation over repeated indices is implied. 

This scalar potential determines the structure of the vacuum, inducing vevs: 

/ rrO\ / ttO\ I -ex "^Rs / - \ ^2 

(i^.) = ^, = = (^.> = ^ (9) 

In particular, the vevs for the right-handed sneutrinos generate effective bilinear couplings: 

h^u LivlHu — Xgv'sHdHu h'jf Li-^Hu — \s-^HdHu = €iLiHu — fiHdHu (10) 
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Since by electroweak symmetry breaking an effective jj, term is generated, it is at the electroweak scale. 
Minimizing the scalar potential gives the following tadpole equations at tree-level 

-l^RsMi^sX; + h.c.) - ^v,{X*XvnsVm + h.c.) - ^vlv,{Xlhl' + h.c.) 

-^v^vns{T^ + h.c.)^0 (11) 

= ^^(.9^ + g'^Wvu + rriH^Vu + ^VuXsX^VBsVm + ^vjvuXsX^ - ^v%^VdiKsXl + h.c.) 
+ ^v,vl^in:hl^ + h.c.) ~ ^VdU^u,iX:hl^ + h.c.) + ^v^v,v,hl%hir 

+ lv^hl%hiyvnsvm - ^VdVRsin + h.c.) + -^v,vr,{TII + h.c.) = (12) 

1^ = ^(5' + g'^)n'v, + i(mi,^. + mi^,)«j - ^Vdvl{X*h}^' + h.c.) 

+ \vr,Vu{hIK-^ + h.c.) - ^VdiXlvRsVRthi* + h.c.) + \vj{vRsVRthi'' (hi*)* + h.c.) 
o 4 4 

+^viv,ihi^{hir + h.c.) + -L=v^vRAn: + ^.c.) = o (13) 

9V 2 1 / X* L ^ 1 * 3 

^ = mj:.a^^VRs - -:VdVuVRs(KsXs + h.c.) + -KsK^Vfi^ 

ClVRs 4 4 

+ ^VuVRsV,{Klhi' + h.c.) + \{vl + vl){XsX;vRt + h.c.) + ^vl[hmhi')*VRt + h.c] 
+ ^v„^Vn[{h'^yh];;'vRt + h.c] - ^VdVj{X;hi'vRt + X*VRthi'^ + h.c.) 
'^^dVuin + h.c.) + -^VuV.iT^^ + h.c.) + J^vj,,vj,^{Tf^ + h.c.) = (14) 



" = vi-vi + vi + vi+vi (15) 



with 



and there is no sum over the index s in Equation (I14p . 

As usual in R-parity breaking models with right-handed neutrinos, see for example the model proposed 
in [s^l, it is possible to explain the smallness of the Vi in terms of the smallness of the Yukawa couplings 
/ly, that generate Dirac masses for the neutrinos. This can be easily seen from Equation (fT3|l . where both 
quantities are proportional. Moreover, as shown in [sH, taking the limit /ij/ — > and, consequently, Vi 0, 
one recovers the tadpole equations of the NMSSM, ensuring the existence of solutions to this set of equations. 



D. Masses of the neutral scalars and pseudoscalars 



In this subsection we work out the main features of the neutral scalar sector mainly focusing on singlets. 
The complete mass matrices are given in Appendix We start with the one generation case which closely 
resembles the NMSSM, considered, for example, in |6l|, [6^. This already implies an upper bound on the 
lightest doublet Higgs mass m{h°), where we will focus on at the end of this subsection. A correct description 
of neutrino physics implies small values for the vevs Vi of the left sneutrinos and small Yukawa couplings hi, 
as we will see later. Neglecting mixing terms proportional to these quantities, the (6 x 6) mass matrix of the 
pseudoscalars in the basis Im{H^, v'^^Vi) given in Appendix El Equation (|A20[) . can be decomposed in 
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two (3 X 3) blocks. By using the tadpole equations we obtain 



MfjH Mis 
Mlo - I (A^^)^ Ml, I (16) 
Af|^, 



with 



Mis = m +^2)"-^- 3^^3, (M|i)„. = i + i (-1),, + [i (5^ + 5'^) , (17) 



where is defined in Equation psp . The parameters f2i are defined as: 

ni^^iXn* + x*k)vI, n2^^{n + n)vR, n,^^{T. + T:)vH (18) 

The upper (3 x 3) block contains the mass terms for Im{Hd), Im[Hu) and I'm(y'^) and we get analytic 
expressions for the eigenvalues: 



m2(P{') = 

2 V^'" '"d v% ) 2 

with r=fi(i7, + f72)f^ + ^ + ^)-V' 



3in^ + n,)nJ^ + ^)-9n^nJ^ + ^) (i9) 



The first eigenvalue corresponds to the Goldstone boson due to spontaneous symmetry breaking. To get 
only positive eigenvalues for the physical states, the condition 

has to be fulfilled, implying that has in general the opposite sign of vr. Additional constraints on the 
parameters are obtained from the positiveness of the squared masses of the neutral scalars. Taking the scalar 
mass matrix from Appendix lAl Equation (jAlip . in the basis Re{H^^H^,v'^ ,Vi) in the same limit as above 
we obtain 

/ Ml^ Mis \ 
= {Mls)^ Mis (21) 
V MlJ 

with 

/\,f2 i\---L i / ' ~ ~^ Vu ~ ~ ~^ ' A'f^ \ ^ ^ ""-^/ 



, - -4 - 



(22) 
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using the additional parameters 

^4 = XX*VdVu > 0, fis = ^KK*vji > 0, = ^ (.9^ + g'^) VdVu > . (23) 

An analytic determination of the eigenvalues is possible but not very illuminating. However, one can use the 
following theorem: A symmetric matrix is positive definite, if all eigenvalues are positive and this is equal to 
the positiveness of all principal minors (Sylvester criterion). This results in the following three conditions 

< {n^ + n2) — + ne— 

Vd Vu 

< (fii + f^a) (^^6 + ^) " + 2f^4^ + 2174f26 - 

< r!3-/2(r!2) , (24) 
where 72(^^2) is given by: 
^1 

/2 (^^2) = with 
^2 



Si = (r^i + ^2) r^s (-2O4 + 21^6) + (o| - 2^4^16) ^5 



Wa ^dJ W Vd 

{ni + ^2) ns^e (^ + ^]+2 (ni + ^2-^4 + 2ne) nl-^ 

2 {2Vli + ^2) (2f^i + 2VI2 - O4 + 21^6) ^4\— + - 

\Vu Vd, 

[l&nl + 8 (4^2 - O4 + Sle) + lOr^ifla (2^22 - + ^e) 

+ ^2 {2^2 - ^4) {2n2 -^4 + 2176)] ^ 



S2 = (r!i + r!2) rie M + -| + 2 (r!i + ^2) {^4 - ^e) 

+ 2n4ne - nl (25) 

The first two conditions are in general fulfilled, but for special values of tan /? or A. Putting all the above 
together we get the following conditions: 

/2(r!2)<r!3</i(r!2) (26) 

It turns out that by taking a negative value of fl^ (oc T^) near f2i^2) one obtains a very light singlet scalar, 
whereas for a value of fis near /i(r22) one gets a very light singlet pseudoscalar. In between one finds a value 
of n^, where both particles have the same mass. This discussion is comparable to formula (37) in ^62^ for 
the NMSSM. Moreover, a small mass of the singlet scalar and/or pseudoscalar comes always together with 
a small mass of the singlet fermion. 

In the n generation case similar result holds as long as and m?c do not have off-diagonal entries 
compared to k. Inspecting Equations (|A15P and (|A24[) it is possible to show that the singlet scalars and 
pseudoscalars can be heavy by appropriately chosen values for the off-diagonal entries of while keeping 
at the same time the singlet fermions relatively light, as will be discussed later. As pointed out in [53], the 
NMSSM upper bound on the lightest doublet Higgs mass of about ^ 150 GeV, which also applies in the 
/ij'SSM, can be relaxed to 0(300) GeV, if one does not require perturbativity up to the GUT scale. 



E. Neutrino masses 



In the basis 



(v^°)^ = (B^w/3^i^,^H°,^.^..0 



(27) 
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the mass matrix of the neutral fermions, see Appendix lA 41 has the structure 

^■■ = (™?o) ■ M 

Here Mh is the submatrix including the heavy states, which consists of the usual four neutralinos of the 
MSSM and n generations of right-handed neutrinos. The matrix m mixes the heavy states with the left- 
handed neutrinos and contains the R-parity breaking parameters. 
The matrix A^„ can be diagonalized in the standard way: 

Mn - U*MnU-^ (29) 

As it is well known, the smallness of neutrino masses allows to find the effective neutrino mass matrix in a 
seesaw approximation 

m^^ = -m^ • MH^'m = -e • m , (30) 

where the matrix ^ contains the small expansion parameters which characterize the mixing between the 
neutrino sector and the heavy states. 

Since the superpotential explicitly breaks lepton number, at least one mass for the left-handed neutrinos 
is generated at tree-level. In the case of the 1 i/'^-model the other neutrino masses are generated at loop- 
level. With more than one generation of right-handed neutrinos additional neutrino masses are generated at 
tree-level, resulting in different possibilities to fit the neutrino oscillation data, see the discussion below. 

1. One generation of right-handed neutrinos 
With only one generation of right-handed neutrinos the matrix ^ is given by 

= KiM - -e,S,3 , (31) 



where the and Ai parameters are defined as 



e^ = -^KvR (32) 
Ai = /iWj -I- etVd (33) 



and Kj^ as 



1 2g'M2^i 



-a 



ith 



" -{vdVuX'' + MRfi) (36) 



Figure 1: Example of one 1-loop correction to the effective neutrino mass matrix involving tlie singlet 
scalar /pseudoscalar. 




Figure 2: 1-loop mixing between gauginos and the right-handed neutrinos. 



and Det{MH) is the determinant of the (5 x 5) mass matrix of the heavy states 

Det{MH) = Im^iX^v^ + AMniiVdVu) - MiAhfiivdVuX^ + Mrij) . (37) 

o 

Using these expressions the tree-level effective neutrino mass matrix takes the form 

imlZh=a^.^j ■ (38) 

The projective form of this mass matrix implies that only one neutrino gets a tree-level mass, while the other 
two remain massless. Therefore, as in models with bilinear R-parity violation [36l . [stI . [63| 1-loop corrections 
are needed in order to correctly explain the oscillation data, which require at least one additional massive 
neutrino. The absolute scale of neutrino mass constrains the A and e parameters, which have to be small. 
For typical SUSY masses order 0(100 GeV), one finds \A\/n^ - lO'^-lO"^ and - lO'^-lO'^. This 
implies a ratio of |e|^/|A| - 10"^-10"^ 

General formulas for the 1-loop contributions can be found in 36] and adjusted to the /xi'SSM with 
appropriate changes in the index ranges for neutralinos and scalars. Important contributions to the neutrino 
mass matrix are due to 6 — 5 and r — f loops as in the models with b-^p [37,] . In addition there are two new 
important contributions: (i) loops containing the singlet scalar and singlet pseudoscalar shown in Figure 
[TJ As shown in [g^, [1^, [la] , the sum of both contributions is proportional to the squared mass difference 
Ai2 = m|j. — mj oc K^vj^ between the singlet scalar and pseudoscalar mass eigenstates. Note that this 
splitting can be much larger than the corresponding ones for the left sneutrinos. Thus the sum of both loops 
can be more important than b — b and r — f loops in the current model, (ii) At loop-level a direct mixing 
between the right-handed neutrinos and the gauginos is possible which is zero at tree-level, see Figure [H 
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2. n generations of right-handed neutrinos 



In this class of models with n > I one can explain the neutrino data using the tree-level neutrino mass 
matrix only. In general one finds that the loop corrections are small if the conditions at the end of this 
section are fulfilled. 

For the sake of simplicity, let us consider two generations of right-handed neutrinos which contains all 
relevant features. The matrix ^ in Equation ([50]) takes the form 



= KiA, + Kia, - -6 



J3 



with 



1 



lis 



The Ka and Ka coefficients are: 



V2 

Aj = ^Vi + eiVd 



1 2g'M2^l 
= —Z. « 



(39) 



(40) 

(41) 
(42) 



Kl 



2gMifi 



8^iDet{MH) 



[vdv'^iMRiXl + Mr2\1) + 2vuMRiMB,2tA 



^5 MR2Xim^ 2 2n 
K^ = —^{Vu - Vd) 



AV2Det{MH] 
MRiX2m^ 



A'2 
Kt 



2gMifi^ 



b 



(m-yV^Vu - AMiM2HVd) 



-V2X2C- 



AV2Det{MH 
The effective neutrino mass matrix reads as 

i'^tZ)v ^ aAiAj + b{A^aj + AjUi) -I- caiUj 



%/2Ai( 



'iDetoVR2 ^ 



(43) 



with 



8V2nDet{MH 
1 



{MRjXlvuVd + MR2X\vuVd + MRjMR,2ti) 

■{""l - Vd){MRlVRiX2 ~ Mr2VR2Xi) 



using Mrs 



-^KsVrs and the determinant of the (6 x 6) mass matrix of the heavy states is 
DetiMn) = ^ [(M^A? + MRiXl){my - SMiNhtiVuVd) + 8MRiMR2Deto] 



with Deto being the determinant of the usual MSSM neutralino mass matrix 

1 2 
—rn.-ji.u,ii].. — A/fi Mnr 



Detci = —m^iiVdVu — M1M2I1 



(44) 

(45) 
(46) 

(47) 

(48) 
(49) 



The mass matrix in Equation (|44p has two nonzero eigenvalues and therefore the loop corrections are not 
needed to explain the experimental data. Two different options arise: 
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• A generates the atmospheric mass scale, a the solar mass scale 

• a generates the atmospheric mass scale, A the solar mass scale 

In both cases one obtains in general a hierarchical spectrum. A strong fine-tuning would be necessary to 
generate an inverted hierarchy which is not stable against small variations of the parameters or radiative 
corretions. Moreover the absolute scale of neutrino mass requires both |A|//i^ and |Q?|//i to be small. For 
typical SUSY masses order 0(100 GeV) we find in the first case |A|//^2 _ iQ-^-iQ-e and \d\/^i - lO-^-lQ-^. 
In the second case we find |A|//i^ ~ 10~®-10~^ and |c?|//J. ~ 10~^-10"^. The ratios including e or a are 
much smaller than those in the 1 v"^ case. We find that 1-loop corrections to ([33]) are negligible if 

<10-3 and C<10-3 (50) 

iAi ^ iAr 

are fulfilled. Note that the mixing of the neutrinos with the higgsinos, given by the third column in the 
matrix ^ in Equation ([39]) . depends not only on a; but also on a. This leads to 1-loop corrections to the 
neutrino mass matrix with pieces proportional to the parameters, as it also happens in the 1 P'^-model. 
Therefore, both conditions in Equation ([50]) need to be fulfilled. Finally, in models with more generations 
of right-handed neutrinos there will be more freedom due to additional contributions to the neutrino mass 
matrix. For example, the case of three generations is discussed in [53J , where the additional freedom is also 
used to generate an inverted hierarchy for the neutrino masses. 



III. CHOICE OF THE PARAMETERS AND EXPERIMENTAL CONSTRAINTS 



In the subsequent sections we work out collider signatures for various scenarios. To facilitate the com- 
parison with existing studies we adopt the following strategy: We take existing study points and augment 
them with the additional model parameters breaking R-parity. These points are SPSla' [67[, SPSS, SPS4, 
SPS9 [111 and the ATLAS SU4 point 0. SPSla' contains a relative hght spectrum so that at LHC a high 
statistic can be achieved, SPS3 has a somewhat heavier spectrum and in addition the lightest neutralino and 
the lighter stau are close in mass which affects also the R-parity violating decays of the lightest neutralino. 
SPS4 is chosen because of the large tan/3 value and SPS9 is an AMSB scenario where not only the lightest 
neutralino but also the lighter chargino has dominant R-parity violating decay modes. In all these points 
the lightest neutralino is so heavy that it can decay via two-body modes, as long as it's not a light v'^. In 
contrast for the SU4 point all two-body decay modes (at tree-level) are kinematically forbidden. As the 
parameters of these points are given at different scales we use the program SPheno 'to'I to evaluate them at 
Q = mz where we add the additional model parameters. Note that we allow /i to depart from their standard 
SPS values to be consistent with the LEP bounds on Higgs masses, discussed below. 

The additional model parameters are subject to theoretical and experimental constraints. In [s^] the 
question of color and charge breaking minimas, perturbativity up to the GUT scale as well as the question of 
tachyonic states for the neutral scalar and pseudoscalars have been investigated . The last issue has already 
been addressed in Section III Dl where we derived conditions on the parameters. By choosing the coupling 
constants A, k < 0.6 in the 1 9'^-model and A^, < 0.5 in the 2 i/'^-model, perturbativity up to the GUT 
scale is guaranteed [52]. Note, that choosing somewhat larger values for A and/or k up to 1 does not change 
any of the results presented below. We also address the question of color and charge breaking minimas 
by choosing A^ > 0, > 0, > 0, T^*" < 0, whereas the Yukawa couplings h]f can either be positive 
or negative, but those values are small < 0(10"^) due to constraints from neutrino physics. Our TJf are 
negative, so the condition (2.8) of [S^l is easy to fulfill. 

Concerning experimental data we take the following constraints into account: 

• We check that the neutrino data are fulfilled within the 2-f7 range given in Table |T] taken from ref. [l2| 
if not stated otherwise. These data can easily be fitted using the effective neutrino mass matrices given 
in Section UTEI 

• Breaking lepton number implies that flavour violating decays of the leptons like /i —^ e"/ are possible, 
where strong experimental bounds exist [3l| . However, in the model under study it turns out that these 
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parameter 


best fit 


2-cr 


Amii[10-5eV2] 
\Aml\[W-^eV^] 

sin^ 012 

sin^ 023 

sin^ 013 


' -"^-0.20 
^■^"-0.11 

m+°"i^ 

"•'-'J^-o.oii 


7.25-8.11 
2.18 - 2.64 
0.27-0.35 
0.39 - 0.63 
< 0.040 



Table I: Best-fit values with 1-a errors and 2-a intervals (1 d.o.f.) taken from |12||. In the following we will refer to 
these angles as 0i2 = 9 sot, 623 = datm and 0i3 = 9r. 



bounds are automatically fulfilled once the constraints from neutrino physics are taken into account 
similar to the case of models with bilinear R-parity breaking [Tlj . 

• Bounds on the masses of the Higgs bosons [3l|, [tI] . For this purpose we have added the dominant 
1-loop correction to the (2,2) entry of the scalar mass matrix in Appendix lA 21 Moreover, we have 
checked in the 1 i>'=-model with the help of the program NMHDECAY ^] that in the NMSSM limit 
the experimental constraints arc fulfilled. 

• Constraints on the chargino and charged slcpton masses given by the PDG [3T| . 

• The bounds on squark and gluino masses from TEVATRON [3l[ are automatically fulfilled by our 
choices of the study points. 

The smallness of the l^p parameters guarantees that the direct production cross sections for the SUSY 
particles are very similar to the corresponding MSSM/NMSSM values. Note that for low values of A the 
singlet states are decoupled from the rest of the particles, leading to low production rates. 

IV. PHENOMENOLOGY OF THE 1 i?'^-MODEL 

In this section we discuss the phenomenology of the 1 P'^-modcl, including mass hierarchies, mixings in the 
scalar and fermionic sectors, decays of the scalar and fermionic states and the correlations between certain 
branching ratios and the neutrino mixing angles. 




20' 

-120 -100 -80 -60 -40 -20 

A, [GeV] 

Figure 3: Masses of the lightest neutralinos x? and the lightest scalar 5*1 = Re{v'^) /pseudoscalar = Im{vl) as a 
function of = T^/k for A = 0.24, k = 0.12, = 150GeV and Ta = 360GeV for SPSla'. The different colors refer 
to the singhno x? (blue), the bino X2 (red), the singlet scalar S? (black) and the singlet pseudoscalar Pf (green). 
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In the following discussion we call a neutralino x'l a bino (singlino) if |A/j+3.ip > 0.5 (|A/j+3,5p > 0.5). As 
discussed below, light scalar or pseudoscalar states appear, especially in case of the singlino being 
the lightest neutralino. In the following we discuss possible mass hierarchies and mixings in more detail. 

The diagonal entry of the singlet right-handed neutrino in the mass matrix of the neutral fermions is 
Mji = -^I'^vr, see Appendix I A 41 A singlino as lightest neutralino is obtained by choosing small values for k 

and/or vr. Since the masses of the four MSSM neutralinos are mainly fixed by the chosen SPS point, we can 
either generate a bino-like or a singlino- like lightest neutralino by varying k and/or vr, where the latter case 
means a variation of A due to a fixed /i-parameter. A light singlet scalar and/or pseudoscalar can be obtained 
by appropriate choices of T\ and T^. An example spectrum is shown in Figure [3) The MSSM parameters 
have been chosen according to SPSla' except for /i = 150 GeV. The scalar state S2 = can easily get too 
light to be consistent with current experimental data, although the production rate e^e^ — > ZS2 is lowered, 
since a mixing with the lighter singlet scalar — v"^ reduces its mass. By reducing ^ the mixing can be 
lowered (see mass matrices) and this problem can be solved. 

Another example spectrum for neutral fermions is shown in Figured] Again SPSla' parameters have been 
chosen, except fx = 170 GeV. As the figure demonstrates for this reduced value of /x the states are usually 
quite mixed, which is important for their decay properties, as discussed below. Note that the abrupt change 
in composition in Xg is due to the level crossing in the mass eigenstates. 




Figure 4: Masses and particle characters of the lightest neutralinos Xi ^ ^ function of k for A — 0.24, fj, — 170GeV, 
Tx = 360GeV and T„ = -k ■ 50GeV for SPSla'. The different colors refer to singlino purity |M+3,5p (blue), bino 
purity lA/;+3,i|^ (red), wino purity [M+aM'^ (black) and higgsino purity |A/i+3,3p + \M+3a\'^ (green). 
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The decay proper ties of the hghtest scalars/pseudoscalars are in general quite similar to those found in 
the NMSSM 0,[6i]. The hghtest doublet Higgs boson similar to the h° decays mainly hke in the MSSM, 
apart from the possible final state 2xi, if kinematically possible. An example is shown in Figure [5l which 
display the branching ratios of 5*2 = h'^ versus m{xi)- Xi in this plot is mainly a singlino (see Figured]), 
variation of k varies its mass, since is kept fixed here. In contrast to the NMSSM this does not lead to 
an invisible Higgs, since the neutralinos themselves decay. For the range of parameters where the decay to 
is large, Xi decays mainly to i/bb, leading to the final state 4 6-jets plus missing energy. Note that the 
Si which is mainly singlet here decays dominantly to bb final states, followed by rr final states. 




Figure 5: Branching ratios Br{S2 ~ hP) as a function of m(xi) for the parameter set of Figure |4] (variation of k). 
The colors indicate the different final states: XiXi (red), bb (blue), t'^t~ (green), cc (orange) and Wqq (brown). 



A. Decays of a gaugino-like lightest neutralino 

We first consider the case of a bino as lightest neutralino. Although to(xi) > fnw in the SPS points we have 
chosen, two-body decay modes are not necessarily dominant. The three-body decay Xi hljV dominated 
by a virtual f also can have a sizeable branching ratio, see Table [ill and Figure [71 The importance of this 
final state can be understood from the Feynman graph shown in Figure [HI giving the dominant contribution 
due to iJ^-Zi-mixing [li — e,/i). 




Figure 6: Dominant Feynman graph for the decay Xi — * Utu with U = e, p. 



In the case li = t there's an additional contribution due to H^-v-mhdng. As Figure [3 shows there 
exist parameter combinations in the A-K-plane, where the decay mode x^ — > liljV is more important than 
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Xi Wl. The strong variation in the branching ratios for SPSla' is mainly due to the strong dependence 
of the partial decay width of Xi ~^ IdjV, where the decays with i — j and i ^ j both play a role. Other 
important final states are Xi ^ Zi/ and in case of a light scalar with to(xi) > m{h'^) the decay Xi ~^ hPv^ 
as demonstrated in Table [TTl 



Brixl) 


SPSla' 


SPS3 


SPS4 


Wl 


23-80 


12-55 


68-72 




11-75 


2-31 


2.6 - 3.9 


Zv 


2.2 - 8.9 


5-28 


25 - 28 


h% 




15-53 


< 2.0 


Decay length [mm] 


1.6-7.0 


0.1-0.5 


1.4 - 1.6 



Table II: Branching ratios (in %) and total decay length in mm of the decay of the lightest bino-like neutralino for 
different values of A G [0.02,0.5] and k G [0.1,0.6] with a dependence of allowed k(A) similar to |53] and to Figure [7] 
and Ta = A ■ 1.5TeV and = ■ lOOGeV. 



^ 0.6 
0.5 



0.4 



0.3 



0.2 



0.1 

0.05 0.10 0.15 0.20 0.25 ^ 

Figure 7: Dependence of allowed k(A) for values of A € [0.02,0.5] and k G [0.1,0.6] and Br{xi Uljv) as function 
of A and k exemplary for SPSla' with ^ = 390GeV, Ta = A • 1.5TeV and = ■ lOOGeV. 

In the /ii^SSM one finds correlation between the decays of the lightest neutralino and the neutrino mixing 
angles, because neutralino couplings depend on the same I/ip parameters as the neutrino masses. Figure [8] 
shows the correlation between the branching ratios of the decay Xi Wl as a function of the atmospheric 
angle. Although a clear correlation is visible it is not as pronounced as in the n generation case, see below 
and [5^, due to inclusion of 1-loop effects in the neutrino masses and mixing angles. 

Also the three-body decay Xi ~^ liljV exemplifies a correlation with neutrino physics. However, this 
decay is connected to the solar angle, see Figure [HI There are two main contributions to this final state: 
Xi Wl — + liljV and Xi ^ t*1 — + liljV. While the former is mainly sensitive to A^, the latter is dominated 
by Ci-type couplings (see Figure [6]), causing the connection to solar neutrino angle. In case the W is on-shell 
as in the SPSla' point, one could in principle devise kinematical cuts reducing this contribution. Such a cut 
can significantly improve the quality of the correlation. 

The SU4 scenario of the ATLAS collaboration [6^ has a very light SUSY spectrum close to the Tevatron 
bound with a bino-like neutralino m(xi) ~ 60 GeV. Thus, for SU4 the lightest neutralino has only three- 
body decay modes. Most important branching ratios are shown in Figure [TUl The lightness of the bino-like 
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10"' 10° IC 

tan^ Oatm 



Figure 8: Ratio versus t&n^ Oatm for different SPS scenarios (SPSla' (black), SPSS (red), SPS4 (green)) 

and for different values of A G [0.02,0.5] and k G [0.1,0.6] with a dependence of allowed k(A) similar to [s^] and to 
Figure [3 and Ta = A • 1.5TeV and T^ = -k- lOOGeV. 
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Figure 9: Ratio versus tan^ 6^80; with same set of parameters as Figure [8] Bino purity |A/'4i|^ > 0.97. To 

the left (a) two-body plus three-body contributions, to the right (b) three-body contributions only. For a discussion 



see text. 



neutralino Xi i^i this scenario implies a larger average decay length of (8 — 90) cm, depending on the parameter 
point in the A-K-plane. Note that the decay length becomes smaller for smaller values of A, n. In general 
the decay length scales as L cx m~^(xi) for m(xi) < mw- Also for this point a correlation between the 
branching ratios and the neutrino mixing angles is found as illustrated in Figure [TT] 

In addition to the SUGRA scenarios discussed up to now we have also studied SPS9, which is a typical 
AMSB point. The most important difference between this point and the previously discussed cases is the 
near degeneracy between lightest neutralino and lightest chargino. This near degeneracy is the reason that 
the chargino decay is dominated by l^p final states. Varying A and k as before we find a total decay length 
of (0.12 - 0.16)mm with Br{xf ^ Wi^) = (42 - 57)%, Br{xf Zl) = (20 - 26)% and Br{x^ h°l) = 
(17 — 40)%. This is especially interesting since, similar to Wl in case of the gaugino-like lightest neutralino, 
the decay to Zl of the chargino is linked to the atmospheric angle, see Figure [T2l 
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Figure 10: Decay branching ratios for bino-like lightest neutralino as a function of k for A G [0.02, 0.5], Tx — \ - 1.5 
TeV, T„ = ■ 100 GeV and for MSSM parameters defined by the study point SU4 of the ATLAS collaboration 
The colors indicate the different final states: liljV (red), Iqiq^ (black), uqq (blue) and (orange). 
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Figure 11: To the left (a) ratio „" — ^""I^' versus tan^ 
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for the SU4 scenario of the ATLAS collaboration 169 



and to the right (b) ratio ^^f:h^'''^'^\ versus tan^ Osoi with same set of parameters as (a). Bino purity |A/'4i|^ > 0.94. 




Figure 12: Ratio ^''^^i ^^'"^ versus tan^ Oatm for the AMSB scenario SPS9 and for different values of A G [0.02, 0.51, 
K G [0.1, 0.6], Ta = A • 1.5 TeV and r„ = -k ■ 100 GeV. 
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B. Decays of a singlino-like lightest neutralino 

We now turn to the case of a singlino-like LSP. As already explained, this scenario is connected to a light 
singlet scalar and pseudoscalar. Recall, that the particles in the fermionic sector are mixed for A, k = 0(10"^) 
due to the reduced /i-parameter as can be seen in Figure SI We will first discuss the average decay length 
of the lightest neutralino Xi- Figure [13] shows the average decay length in meter for different SPS scenarios 
as a function of the mass of the lightest neutralino m{xi)- Composition of the neutralino is indicated by 
colour code, as given in the caption. A, k, and are varied in this plot. Note that by variation of T„ the 
parameter points in Figure [T3] are chosen in such a way, that all scalar and pseudoscalar states are heavier 
than the lightest neutralino. Singlino purity in this plot increases with decreasing mass and for pure singlinos 
the decay length is mainly determined by its mass and the experimentally determined neutrino masses. For 
neutralino masses below about 50 GeV decay lengths become larger than 1 meter, implying that a large 
fraction of neutralinos will decay outside typical collider detectors. Note that if one allows for lighter scalar 
states so that at least one of the decays Xi ~^ ^^{Pi)^ appears, the average decay length can be easily 
reduced by several orders of magnitude. 




20 40 60 80 100 120 140 

?Ti(x;) [GeV] 

Figure 13: Decay length of the lightest neutralino Xi in m as a function of its mass m(xi) in GeV for different 
values of A G [0.2,0.5], k G [0.025,0.2] and fi G [110, 170]GeV with a dependence of allowed k(A) similar to ^ 
and to Figure [7] and Tx — \ ■ 1.5TeV, whereas G [—20, — 0.05]GeV is chosen in such a way, that no lighter scalar 
or pseudoscalar states with {m{Si),m{Pi)} < m(xi) appear. Note that the different colors stand for SPSla' (real 
singlino, |A/'45|^ > 0.5) (gray), SPSla' (mixture state) (black), SPSS (real singlino) (blue), SPSS (mixture state) (red) 
and SPS4 (mixture state) (green). 

Again typical decays are Wl, Iqiqj, Zv, vqq, liljV and the invisible decay to iv. For the region of m{x\) 
below the W threshold see Figure fT4l The dominance of vbb for smaller values of m(xi) is due to the decay 
chain x\ — > S^v — > vhb, whereas for larger values of 'ti(xi) we find m{Si) > m{xi)- Final state ratios show 
correlations with neutrino physics also in this case. As an example we show liljh' branching ratios versus 
the solar neutrino mixing angle in Figure [T5l Singlino purity for this plot |A/'45|^ G [0.75,0.83] and mass 
iTi{xi) S [22, 53] GeV. The absolute values for the branching ratios are comparable to those of the described 
SU4 scenario with a bino-like lightest neutralino. We note that for the parameters in Figure [12] the light 
Higgs = decays to XiXi with a branching ratio of Br{S^ = /i° -> XiX?) = (21 - 91)%. 

Up to now we have considered values of A and k larger than 10^^. For very small values of these couplings, 
the singlet sector, although very light, effectively decouples. This implies that R-parity conserving decays 
of X2j 6.g. decays to final states like Xi'5'i , XiPii Xi^"*"^" or XilQi £^re strongly suppressed and the decay 
modes dominate, implying decays with correlations as in the case of the explicit h-l^p . 
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Figure 14: Singlino decay branching ratios as a function of its mass, for the same parameter choices as in Figure [5] 
The colors indicate the different final states: v\A) (blue), liljV (red), Iqiq^ (black), ?>v (orange) and vqq {q ^ b, green). 
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Figure 15: Ratio gl.[fCl7j) versus tan^ 9soi for the SPSla' scenario and A £ [0.2, 0.5], fi £ [110, 170]GeV, k = 0.035, 
Ta = A • 1.5 TeV and = -0.7 GeV. 



V. PHENOMENOLOGY OF THE n P'^-MODEL 



In the previous section the phenomenology for the one generation case of the model has been worked out 
in detail. Most of the signals discussed so far are independent of the number of right-handed neutrinos. 
However, the n generation variants also offer some additional phenomenology, which we discuss here for the 
simplified case of n = 2. 

In a model with one right-handed neutrino superfield a light singlino will always imply a light 
scalar/pseudoscalar. This connection between the neutral fermion sector and scalar/pseudoscalar sector 
is a well-known property of the NMSSM (see again [6l|, [s^])- In models with more than one generation of 
singlets, the off-diagonal terms in Equation ([5]) induce mixing between the different generations of singlet 
scalars and pseudoscalars. This opens up the possibility, not considered in previous publications 51, 52, 53], 
to have the singlet scalars considerably heavier than the singlet fermions. 

Let us illustrate this feature with a simple example. Imagine a light singlino vf, and a heavy singlino 
1/2, in a model with non-zero trilinear couplings T^^^. In that case, the contributions to the mass of the 
vf, scalar or pseudoscalar, coming from the large value of vii2 are proportional to T^^^. Without these 
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Figure 16: Masses of the scalar states Re{ui) (green), Re{i>2) (red) and ft" (blue) and the pseudoscalar states Im(ui) 
(dashed green), Jm(i>f ) (dashed red) and Im{i>i) (dashed blue) as a function of vr2 for different values of T;^^^ = T^^^. 



To the left (a) r„ ^ = 



whereas to the right (b) 



-2 GeV. The MSSM parameters have been 



taken such that the standard SPSla' point is reproduced. The light singlet parameters rei = 0.16 and vm = 500 GeV 
ensure that in all points the lightest neutralino is mostly uf, with a mass of 47 - 48 GeV. In addition, = 300 GeV 
and Tl G [10, 200] GeV. 



contributions the mass of iyf would only depend on the small v^i, thus making it light like the singiino of 
the same generation. With non-zero T^^^ the mass of both i/f are dominated by the larger of the vjis- This 
feature is demonstrated in Figure [TBI In the two plots the lightest neutralino is mostly with a mass of 
~ 50 GeV. These plots show the dependence of the masses of the singlet scalar states Re{Df) and i?e(i>2) 
and the corresponding pseudoscalar states Im{i>f) and Im(i>2) with vji2 for different values of T^^^ = T^^^. 
The masses of the light Higgs boson and the lightest left-handed sneutrino /to(i>i) are also shown for 
reference. Note that for T^^^ = T^^^ = the mass of the state 7?e(£'J) does not depend on vm, whereas for 
T^^^ = T^^^ = —2 GeV the lightest singlet scalar becomes heavier for larger values of vr2- The same feature 
is present in the pseudoscalar sector, where the effect is even more pronounced. 



A. Correlations with neutrino mixing angles in the n z/'^-model 

The connection between decays and neutrino angles is not a particular property of the 1 P'^-model and is 
also present in a general n 9'^-model. However, since the structure of the approximate couplings x? — — If 
is different, see Appendix [B] we encounter additional features for n = 2. 

As explained in Section HIE 21 we have now two possibilities to fit neutrino data. If the dominant con- 
tribution to the neutrino mass matrix comes from the A^Aj term in Equation (|44[) one can link it to the 
atmospheric mass scale, using the aiaj term to fit the solar mass scale. This case will be called option fitl. 
On the other hand, if the dominant contribution is given by the aiaj term one has the opposite situation, 
where the atmospheric scale is fitted by the parameters and the solar scale is fitted by the A; parameters. 
This case will be called option fit2. 

For the case of a bino-like lightest neutralino one can show that the coupling is proportional to A^ 
whereas for the case of a singiino- like lightest neutralino the dependence is on a^, as shown in Appendix 
m Figure [17] shows the ratio Br{xi W^j,)/Br{x^ Wt) versus tun^ieatni) (left) and Br{xi 
We)/ Br{x1 W^)"^ + Br{x1 WtY versus sin^(0fl;) (right) for a bino LSP and option fitl. The correla- 
tion with the atmospheric angle and the upper bound on Br{x1 — > We)/ Br{xi ~^ Wiu,)"^ + Br{x'( WtY 
from SIT? [Oil) is more pronounced than in the 1 j)'^-model, because we fit neutrino data with tree- level physics 
only. Recall that this implies that the ratio |e|^/|A| is much smaller than in the plots shown in the previous 
section. A correlation between Br{xi We) / ^ Br{xi — > W^)^ + Br{xi Wt)^ and tan^ (6'so;) is found 
instead, if neutrino data is fitted with option fit2, as Figure [T51 shows. 
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Figure 17: To the left (a) ratio '^^[^Zwr] versus tan^{ea 

versus sin^(&j;) for a bino LSP. Bino purity 1^^411^ > 0.9. Neutrino data is fitted using option fitl. 
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Figure 18: Ratio Br^xl-'W.) 

^Br{x\^Wti.f+Br(x<l^WTy 

data is fitted using option fit2. 



versus tan^(&ao;) for a bino LSP. Bino purity \Nii\^ > 0.9. Neutrino 




Figure 19: Ratio „ Br(x';^w.) 

y/Br{x^l^Wi^)2+Br(x'>^WT)2 

Neutrino data is fitted using option fitl. 
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versus tan'^(0ao;) for a singlino LSP. Singlino purity |A/'45p > 0.9. 
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Figure 20: To the left (a) ratio g^j^o^^^j versus tan'^ (Satm) and to the right (b) ratio 
versus sin^(&_R) for a singlino LSP. Singlino purity \Af45 
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^ > 0.9. Neutrino data is fitted using option fit2. 



For the case of a singlino LSP the correlations and types of fit to neutrino data are swapped with respect 
to the gaugino case. Since the couplings Xi ~ — If are mainly proportional to a^, instead of A^, a scenario 
with a singlino LSP and option fitl (fit2) will be similar to bino LSP and option fit2 (fitl). This similarity is 
demonstrated in Figures [TOl and EOl To decide which case is realized in nature, one would need to determine 
the particle character of the lightest neutralino. This might be difficult at the LHC, but could be determined 
by a cross section measurement at the ILC. We want to note, that in the 2 P'^-model we cannot reproduce 
all correlations for a singlino LSP presented for the 3 P'^-model in [ssl ]. 

The results shown so far in this section were all calculated for the SPSla' scenario. We have checked 
explicitly that for all the other standard points results remain unchanged. We have also checked that for 
a LSP with a mass below raw the three-body decays Xi ~^ ^'QiQjj mediated by virtual W bosons, show the 
same correlations. 

A final comment is in order. In a n 9'^-model with n > 2, the effective neutrino mass matrix will 
have additional terms with respect to ()44p . due to the contributions coming from the new right-handed 
neutrinos. For this richer structure there is one additional contribution to which could be sub- 

dominant. Therefore, one can imagine a scenario in which a third generation of singlets produces a negligible 
contribution to neutrino masses while the corresponding singlino, v^, is the LSP. In such a scenario the 
correlations between the LSP decays and the neutrino mixing angles will be lost. 



B. xi decay length and type of fit 

As already discussed we have two different possiblities to fit neutrino data: A generates the atmospheric 
mass scale and a the solar mass scale (case fitl), or vice versa (case fit2). It turns out that the decay length 
of the lightest neutralino is sensitive to the type of fit, due to the proportionality between its couplings with 
gauge bosons and the l^p parameters (see Appendix [B] for exact and approximated formulas of the couplings 
Xi — — If and their simplified expressions in particular limits). For example, a singlino- like neutralino 
couples to the gauge bosons proportionally to the ai parameters. This implies that its decay length will 
follow L (X and obeys the approximate relation 

L{fit2) ~ m,oi ~ ' ^ ' 

In Figure!^ the decay length of the lightest neutralino and its dependence on the type of fit to neutrino data 
is shown. Once mass and length are known this dependence can be used to determine which parameters 
generate which mass scale. Note that this feature is essentially independent of the MSSM parameters. 
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Figure 21: Decay length of the lightest neutralino and its dependence on the type of fit to neutrino data. To the 
left (a) the decay length of the lightest neutralino versus for the case fitl (red) and the case fit2 (blue). To 

the right (b) the ratio L(fitl)/I/(fit2) versus m(x?). The MSSM parameters have been taken such that the standard 
SPSla' point is reproduced. The light singlet parameter k is varied in the range k € [0.01, 0.1]. In all the points the 
lightest neutralino has a singlino purity higher than 0.99. 



However, this property is lost if either the lightest neutralino has a sizeable gaugino/higgsino component or 
if there are singlet scalars/pseudoscalars lighter than the singlino. 



C. Several light singlets 

In scenarios with two (or more) light singlets, the phenomenology has additional features. The light 
Higgs boson /i" can decay with measurable branching ratios to pairs of right-handed neutrinos of different 
generations. Similarly, the bino can decay to the different light right-handed neutrinos. 

In the following, the case of two light singlinos and two light scalars/pseudoscalars will be considered. For 
the neutral fermion sector this implies that the mass eigenstates Xi and X2 will always be the singlets i^f 
and a-nd the bino will be the Xs- In the scalar sector one has two very light mostly singlet states 5'° and 
S2, which are consistent with the LEP bounds. Finally, the state will be the light doublet Higgs boson 
One can also have light singlet pseudoscalars. 

The decays of a bino-like Xg can be very important to distinguish between the one generation model and 
models with more than one generation of singlets. In principle, the most important decay channels strongly 
depend on the couplings of the bino to the two generations of singlinos and the configuration of masses 
of singlinos and scalars. Therefore, a general list of signals cannot be given. Nevertheless, there are some 
features which are always present: 

When kinematically allowed, the decays X3 ^ Xi 2 ^li^i) dominate, with the sum of the branching ratios 
typically larger than 50 %. The relative importance of the different channels is mainly dictated by kinematics. 
This feature is illustrated in Figure [HI where these two quantities are shown as a function of the mass of 
the lighest neutralino. The MSSM parameters are fixed to the standard point SPSla', with light singlet 
parameters taken randomly. One can see that the relative importance of each singlino cannot be predicted 
in general, but both branching ratios are at least of order 10^"^ — 10^"*, given enough statistics. For very light 
singlinos two-body decays including scalars and pseudoscalars are open, and thus both Br{x3 —> Xi) and 
Br{x^ X2) are close to 50%, as expected if the values of the singlet parameters are of the same order for 
the two light generations. On the other hand, if the mass of the lightest neutralino is increased some of the 
two-body decays are kinematically forbidden, specially those of the X2j which has to be produced through 
three-body decays, leading to a suppresion in Br{x'^ x^)- Note that it is also possible to find points where 
the decay mode X3 — > Xi 2 '^2(^2) ^las a branching ratio about 10%-20%, giving additional information. 

The other possible signals are the usual bino decays of the NMSSM. Final states with standard model 
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Figure 22: Branching ratios Br{x3 ~ o Xi) (red) and Br{x3 = o ~* X2) (blue) as a function of the mass 
of the lightest neutralino for the scenario considered in Section IV CI The MSSM parameters have been taken such 
that the standard SPSla' point is reproduced, whereas the singlet parameters are chosen randomly in the ranges 
vm,VR2 G [400,600] GeV, Ai,A2 G [0.0,0.4], T^" = T^^z e [-15,-1] GeV, T^'^ = T^^^ G [-1.5,-0.005] GeV and 
T\,T^ G [0, 600] GeV. k,\ = K2 — 0.16 is fixed to ensure the lightness of the two singlinos. 
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Figure 23: Higgs boson decays as a function of the mass of the lightest neutralino for the scenario considered in 
Section rV CI To the left (a) the standard decay channel hP — > bb, whereas to the right (b) the exotic decays to pairs 
of singlinos XiXi (red), X1X2 (blue) and X2X2 (black). The parameters are chosen as in Figure [52] 



particles, like x\ 2^^^^ Xi 299i become very important when the decays to scalars and pseudoscalars are 
kinematically forbidden. 

In addition, the decays of the light Higgs boson hP can also play a very important role in the study of 
the different generations, provided it can decay to final states including Xi or X2- In this case typically the 
standard Higgs boson decays are reduced to less than 40%, completely changing the usual search strategies. 

In Figure [53] the branching ratios of standard and exotic Higgs boson decay channels are shown. The left 
plot shows the suppressed branching ratio of the standard 66 channel. The main decay channel is Xi Xi: but 
there is a sizeable branching ratio to x5 X2- Note that X2 decays dominantly to x5 plus two SM fermions. 
This feature allows us to distinguish between the 1 P'^-model and models with more than one generation of 
singlets. Finally, the branching ratio to X2 X2 small due to kinematics, but leads to interesting final states 
with up to eight b-jets plus missing energy. 

A final comment is in order. In these kind of scenarios with many light singlets Xi decays to z^66 can be 
dominant. This will reduce the available statistics in the interesting liljV and Iqiq^ channels. Moreover, the 
correlations are less pronounced due to mixing effects in the singlet sector. 
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VI. DISCUSSION AND CONCLUSIONS 

We have studied the phenomenology of the /zz/SSM. This proposal solves at the same time the /z-problem of 
the MSSM and generates small ncntrino masses, consistent with data from neutrino oscillation experiments. 
Neutrino data put very stringent constraints on the parameter space of the model. Both the left-sneutrino 
vacuum expectation values and the effective bilinear parameters have to be small compared to MSSM soft 
SUSY breaking parameters. As a result all SUSY production cross sections and all decay chains are very 
similar to the NMSSM, the only, but phenomenologically very important, exceptions being the decay of the 
LSP and NLSP (the latter only in some parts of the parameter space) plus the decays of the lightest Higgses. 

We have discussed in some details two variants of the model. In the simplest version with only one gener- 
ation of singlets 1-loop corrections to the neutralino-neutrino mass matrix need to be carefully calculated in 
order to explain neutrino data correctly. The advantage of this minimal scheme is that effectively it contains 
only six new (combinations of) .^p parameters, which can be fixed to a large extend by the requirement that 
oscillation data is correctly explained. This feature of the model is very similar to explicit bilinear R-parity 
breaking, although, as we have discussed, the relative importance of the different 1-loop contributions is 
different in the /xz^SSM and in bilinear I^p . Certain ratios of decay branching ratios depend on the same pa- 
rameter combinations as neutrino angles and are therefore predicted from neutrino physics, to a large extend 
independent of NMSSM parameters. We have also calculated the decay length of the LSP, which depends 
mostly on the LSP mass and the (experimentally determined) neutrino masses. Lengths sufficiently large to 
obscirvc displaced vertices are predicted over most parts of the parameter space. However, for neutralinos 
lighter than approximately 30 GeV, decay lengths become larger than 10 meter, making the observation of 
I^p difficult for LHC experiments. However, if there is a singlet scalar or pseudoscalar with a mass smaller 
than the lightest neutralino, Xi — * ^miPm)^ '^'^ dominant decay mode and the corresponding decay 
lengths become much smaller, such that the displaced vertex signature of I^p might even be lost in some 
points of this part of parameter space. On the other hand, in case the mass of the lightest scalar is larger 
than twice the singlino mass, the decay 5"^ 2x? becomes important, both for S'JJj ~ v"^' and ~ hP . If 
this kinematical situation is realized also the Higgs search at the LHC will definitely be affected. 

The more involved n generation variants of the /xi/SSM can explain all neutrino data at tree-level and 
therefore are calculationally simpler. Depending on the nature of the neutralino, neutralino LSP decays 
show different correlations with either solar or atmospheric neutrino angles. This is guaranteed in the two 
generation version of the model and likely, but not always true, for n generations. If the NMSSM coupling A 
is sufficiently small also the NLSP has decays to I^p final states with potentially measurable branching ratios. 
In this part of parameter space it seems possible, in principle, to test both solar and atmospheric neutrino 
angles. If only the singlino(s) are light, i.e. the singlet scalars are heavier than, say, the hP , the decay 
length of the singlino is very sharply predicted as a function of its mass and either the solar or atmospheric 
neutrino mass scale. If both, singlinos and singlet scalars (or pseudoscalars) are light, bino NLSP and hP will 
decay not only to the lightest singlinos/singlets but also to next-to-lightest states. This leads to enhanced 
multiplicities in the final states and the possibility to observe multiple displaced vertices. 

We now briefly discuss possible differences in collider phenomenology of the jivSSM. and other R-parity 
breaking schemes. Different models of R-parity breaking appear clearly distinct at the Lagrangian level. 
However, at accelerator experiments it can be very hard to distinguish the different proposals. This can be 
easily understood from the fact that for a heavy singlet sector all I^p models approach necessarily the MSSM 
with explicit R-parity breaking terms. It is therefore an interesting question to ask, what - if any - kind of 
signals could exist, which at least might hint at which model is the correct description of IjLp . Given the 
large variety of possibilities and the very limited predictive power of the most general cases, any discussion 
before the discovery of SUSY must be rather qualitative. 

First one should mention that not all I]ip models explaining neutrino data show correlations between LSP 
decay branching ratios and neutrino angles. Especially the large number of free parameters in trilinear 
models exclude the possibility to make any definite predictions, .^p models which do show such correlations, 
on the other hand, lead usually to very similar predictions for the corresponding LSP decays. For example, 
fitting the atmospheric data with tree-level I^p terms, a bino LSP in explicit bilinear models and in the 
/Lti^SSM decay with the same ratio of branching ratios into Wl (or Iqiqj) final states. Thus, to distinguish 
the different proposals other signals are needed. 

We will briefly discuss the main differences in collider phenomenology between the following three propos- 
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Displaced vertex 


Comment 


Br(invisible) 


Higgs decays 




Yes 


Visible 


< 10 % 


standard 




Yes/No 


anti-correlates with invisible 


any 


non-standard (invisible) 




Yes/No 


anti-correlates with non-standard Higgs 


< 10 % 


non-standard 



Table III: Comparison of displaced vertex signal, completely invisible final state branching ratios for LSP decays and 
lightest Higgs decays for three different R-parity violating models. For a discussion see text. 

als: (i) MSSM with explicit bilinear terms (b-^p ); (ii) Spontaneous l^p {s-^p ) model and (iii) /ii^SSM. Table 
pill) shows a brief summary of this comparison. Differences occur in (a) the observability of a displaced 
vertex of the lightest neutralino decay; (b) the upper limit on the branching ratio of the lightest neutralino 
decaying completely invisible and (c) standard versus non-standard lightest Higgs decays. 

The decay length of the lightest neutralino is fixed in both, the h-^p model and the ^i^SSM, essentially 
by the mass of the lightest neutralino and the experimentally determined neutrino masses. For m{xi) larger 
than the W-mass decay lengths are typically of the order of 0{mm) and proportional to ?7i~^(xi)- For hghter 
neutralinos, larger decay lengths are expected, see Figures [T3] and [51] which scale like m~^(xi). Shorter 
decay lengths are not possible in h-I^p and possible in the /ii^SSM only if at least one (singlet) scalar or 
pseudoscalar is lighter than x\, when x° S^{P^-i)v dominates. Since in the /ij/SSM the singlet scalars 
decay with a short decay length to 66, one expects that in the /xz^SSM short Xi decay lengths correlate with 
the dominance of bb -f missing energy final states. In the s-^p , on the other hand, the Xj* decay length can 
be shorter than in the b-^p , due to the new final state ji^ ^ J -\- v, where J is the Majoron. Therefore, 
different from the /i^SSM, the neutralino decay length in the s-^p model anti-correlates with the branching 
ratio for the invisible neutralino decay. 

Finally, in the b-^^pOne expects that the decay properties of the lightest Higgs (/i") are equal to the MSSM 
expectations, the only exception being the case when — > 2xi is possible kinematically, in which the Xi 
decays themselves can then lead to a non-standard signal in the Higgs sector. This is different in s-.^p , where 
for a low-scale of spontaneous R-parity breaking, the can decay to two Majorons, i.e. large branching 
ratios of Higgs to invisible particles are possible. In the /^I'SSM the decays can be non-standard, if the 
lightest singlino is lighter than m{h^)/2. However, since the singlinos decay, this will not lead to an invisible 
Higgs, unless the mass of the singlino is so small, that the decays occur outside the detector. 

To summarize this brief discussion, b-^p , s-^p and /i^^SSM can, in principle, be distinguished experi- 
mentally if the singlets are light enough to be observed in case of s-^p and /xz^SSM. We note in passing that 
we have not found any striking differences in collider phenomenology of the /xi/SSM and the NMSSM with 
explicit bilinear terms. 

In conclusion, the /ii^SSM offers a very rich phenomenology. Especially scenarios with light singlets deserve 
further, much more detailed studies. 
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Appendix A: MASS MATRICES 

In the scalar mass matrices shown below the tadpole equations have not yet been used to reduce the 
number of free parameters. 

1. Charged Scalars 



In the basis 



{S-y = iH~,{H+)*,iL,fiL,fL,e*ji,ii*R,f*R) (Al) 
the scalar potential includes the term 

V D (5-')^M|±5+' , (A2) 
where M|± is the (8 x 8) mass matrix of the charged scalars. In the ^ = gauge it can be written as 



Mir M?r 



^5. = I ^ I • (A3) 



The (2 X 2) Mfjjj matrix is given by: 



+ ^KXtVRsVRt + -^Vi{hEh^E)ij'"j 



{M^nh, = (Mi,):, 

(Mi,),, = rnl^ + l[ig^+g'^)vl + i9'-g''){vl+v',+vl+vi)] 

+ -XsKvRsVRt + ^VRsVRthl'ih^)* (A4) 
The (6 X 2) matrix that mixes the charged Higgs bosons with the charged sleptons is 

<r=f50 (A5) 



with: 



M'hr 



{^hhta = \9^VdVi-]^XlK'^VRsVRt-^Vd{hEh}E)ij'"3 

{Mhr)^^ = -\vuVRs{h*E)J^hi' - ^Vj{T^)ji 

{^HR)i2 = -\>^sVRsVj{h*E)ji - ]^Vd{hE)jiKvRs (A6) 
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Finally, the (6 x 6) mass matrix of the charged sleptons can be written as 



with: 



^1= r^f" "^Im (A7) 



(MDij = {ml)^^ + l{9'^-9^Kvl-vl + vl+vi + vl)6,j + ^g\vj 
+ lvl{hEhi).. + ^VRsVRthi'ihi*)* 

+ \vl{h\hE)^. + \vkVm{h^E).^{hE)^. (AS) 



2. Neutral Scaleirs 



In the basis 

[S'^'Y = Re{HlHlvl,yi) (A9) 

the scalar potential includes the term 

V D (50')^Mjo50' (AlO) 
and the ((5 + n) x (5 + n)) neutral scalar mass matrix can be written as 

Mh = {MlsT Ml, . (All) 

V {Ml-f [Mlf Ml- I 

The matrix elements are given as follows: 

(M^i/)n = ^H, + \{9^ + 9'^){^vl -vl+vl+ vl + vl) 

+ ^KKvRsVRt + ^vlXsX* 
(Mhh) 12 = -\{9^ + 9'^)vdVu + XsK'^dVu - ^Vrs{>^sK*s + h.c.) 

-\vuVi{KK + h.c.) - .^VRsin + h.c.) 



{MU)2i = {M%h) 



21 V'-HH) 12 

{M%h)22 = - ^(.9" + g'^M - ^< +vl + vl + vl) 

+ \\sKvRsVRt + ^vl\,Xl + ^VRsVRth}J{h}^)* + ^v^VjiKyhi' 
-\vdVi{KK + h.c.) (A12) 
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The matrix elements are given as follows: 

{mU)^2 = l<si^s< + h-c.) + ^VRsm + h.c.) 

Mhh)^^ = mij^-\{9^ + g'^M-vl + vl + vl + vl) 

+\Xs\*tVRsVm + \vlXsXl + ^VRsVmKi^iT + \viVj{K)*K' 
-]^VdVi{X*h}j + h.c.) (A21) 



{M%s)ts = -\vuVrs{XIks + h.c.) + \'^ViVRt{Xlh^^ - Xthi' + h.c.) 

+ ^^'d(rA^ + h.c.) - ^ViiUl + h-c.) (A22) 



--vl{X*h}: + h.c.) - -VRsVm{X*X + h.c.) 



Ml),, 



-^vkKhl^ + h.c.) ^VRsiTll + h.c.) 



(A23) 



Mss)st = ^[("^iO^t + ("^ic)t«] + \{XsX; + h.c.){vl + vl) + -^VdVu{XlKs + h.c.)5st 
+^KsK*yR,5st - \vuVi{K*,h}j + h.c.)5st + \vl[{K'J)*K-^ + h.c] 

+ \v,vAiKyhi' + h.c] - \vdvAX:h^^ + XtiK)* + h.c] 

--l^VRu{n''' + h.c.) (A24) 



Mis) si = \vuVrs{k*MJ + h.c.) + \'Y^ VdVRtiXthi' - A>t* + h.c.) 

+1 E ^j^ntihi'ihir - hi\hir + h.c] - ^Mn: + h.c.) (A25) 



+\vl[K\hiy + h.c] + \vRsVRt[K\hi'r + h.c] (A26) 



4. Neutral Fermions 



In the basis 
the lagrangian of the model includes the term 

with the ((7 + n) x (7 + n)) mass matrix of the neutral fermions, which can 



Mr mo 



D 







M^o is the usual mass matrix of the neutralinos in the MSSM 



/ Ml 



-\9'vd \g'vu \ 







Mo 



-\g'vd \gvd 
i„'„, 1 



Igvd 





-\gvM 



gvu 



/ 



with 



The mixing between the neutralinos and the singlet is given by 
m^o^, is the neutralino-neutrino mixing part 

■^g'vi ei 

= I -hg''"'2 \gv2 €2 



with 



The neutrino Dirac term is 



and finally Mr is 



-hg'vi \gvz €3 



1 " 



{Mr) St = —7=KaVRs5at ■ 

v2 
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5. Charged Fermions 

In the basis 

{4,+)^ ^(w+,H+,e^,ti^,T^), (A37) 
the (5 X 5) mass matrix of the charged fermions is given by 



/ M2 













\ 




A* 












-ei 














-£2 














-£3 










/ 



Approximate formulas for the coupUng ^ — 
interaction lagrangian 



Appendix B: COUPLING x? -W^ - if 

If can be obtained from the general yf^ 



(Bl) 



where 



Ocnw 
Lil 



Ocnw 
Ril 



y'^Lji ) 



(B2) 



The matrix M diagonalizes the neutral fermion mass matrix (see Appendix (|A 4p ) while the matrices Li and 
V diagonalize the charged fermion mass matrix (see Appendix (|A 5|) ). 

As was already mentioned for the case of neutral fermions in Section IIIEi it is possible to diagonalize the 
mass matrices in very good approximation due to the fact that the I^p parameters are small. Defining the 
matrices ^, and that will be taken as expansion parameters, one gets the leading order expressions 



u = 



V = 



where I3 is the (3 x 3) identity matrix. The expansion matrices and ^^j are 




(B3) 



V2Det+ 



gvdh 



E 



2Det+ 

V^Vdh 
2DeU 



{2i? + g^vl) A, 



E 



Mae, 



2^Det+ 

{vdjJ' + MiVu) Aj 



(B4) 
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where Det+ = — ^g^VdVu + Af2M is the determinant of the MSSM chargino mass matrix, /i = -^X^Vfis and 
ti = -^VRsh]^ ■ The expressions for the matrix ^ depend on the number of singlet generations in the model. 
Particular cases can be found in (|3ip and ([55)1 . 

Using the previous equations and assuming that all parameters are real , one gets the approximate formulas 

^cnw _ 9 r9Nl2M , S^^^W /V A 1 

V2 Det+ /I 2fiDet+ ^ 
Omi = 2^(/^^)»^^L + 2^10^^; M 

Lil — \^Lil ) 

O'i^^i = (OSr)* • (B5) 

It is important to emphasize that all previous formulas, and the following simplified versions, are tree-level 
results. More simplified formulas are possible if the lightest neutralino has a large component in one of the 
gauge eigenstates. These particular limits are of great interest to understand the phenomenology: 



Bino-like Xi 

This limit is caracterized by N^i = 1 and Nim = for m 7^ 1. One gets 

= . (B6) 

For the 1 9'^-model this implies that a bino-like x? couples to WU proportionally to A^, see Equation (|3ip . 
without any dependence on the ti parameters. 

On the other hand, for the 2 P'^-model, the more complicated structure of the ^ matrix, see Equations 
([5^ and implies a coupling of a bino-like x\ with WU dependent on two pieces, one proportional to 
h-i and one proportional to : 

6i = ^^^(aA,-F6a,) (B7) 

However, a simple estimate of the relative importance of these two terms is possible. By assuming that all 
masses are at the same scale msuSY, the couplings k and A are of order 0.1, and the terms h]^ and Vi are 
of order hjjt^ and msusYhj^^ respectively, one can show that aAi 200 foa^. Therefore, one gets a coupling 
which is proportional, in very good approximation, to A^, as confirmed by the exact numerical results shown 
in the main part of the paper. Similar arguments apply for models with more generations of right-handed 
neutrinos. 

In conclusion, for a bino-like neutralino the coupling Xi ~ ~ proportional to A^ to a good 

approximation. 



Higgsino-like Xi 

This limit is caracterized by Nf^ + Nf^ = 1 and Nim = for to 7^ 3, 4. If the coupling O'j^^ is neglected 
due to the supression given by the charged lepton Yukawa couplings, one gets 

^^'^ ' "V2 

^ . (B8) 

Equations (PT|) and ([55)1 show that the terms cancel out in the coupling (jB8[) . and therefore one gets 
dependence only on A^ in the 1 P'^-model, and {Ai,ai) in the 2 ;7'^-model. However, this cancellation is not 
perfect in 0^5" ^^"^ thus one still has some dependence on e^. 
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Singlino-like Xi 

The limit in which the right-handed neutrino ly^ is the hghtest neutraUno is caracterized by Nij^ — 1 for 
TO > 5 and Nn = for / 7^ to. One gets 



9_ 

= . (B9) 

For the 1 9'^-model this expression impUes that a pure singUno-Hke Xi couples to Wh proportional to A^, 
see Equation pip , without any dependence on the parameters. This proportionality to is different to 
what is found in spontaneous R-parity violation, where the different structure of the corresponding 1^ matrix 
[73! impUes that the singlino couples to Wk proportionally to e^. 

For the n P'^-model one finds that the coupling Xi ~ ~ If for £^ singlino-like neutralino has little 
dependence on A^. For example, in the 2 P'^-model one finds that the element ^^5, corresponding to the 
right-handed neutrino uf, is given by 

= T^^>7T777^("« ~ - V 2A2C -1- — 2-6 a, . (BIO 

The coupling has two pieces, one proportional to Ai and one proportional to ai. However, the ai piece gives 
the dominant contribution, as can be shown using an estimate completely analogous to the one done for 
a bino-like Xi- In this case, the ratio between the two terms in Equation (jBlOp is ai-piece ~ 8 A^-piece, 
sufficient to ensure a very good proportionality to the ai parameters. This estimate has been corroborated 
numerically. 
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